Spectral function of electron-phonon models by cluster perturbation theory 
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Cluster perturbation theory in combination with the Lanczos method is used to compute the 
one-electron spectral function of the Holstein polaron in one and two dimensions. It is shown 
that the method allows reliable calculations using relatively small clusters, and at the same time 
significantly reduces finite-size effects. Results are compared with exact data and the relation to 
existing work is discussed. We also use a strong-coupling perturbation theory-equivalent to the 
Hubbard I approximation-to calculate the spectral function of the quarter-filled Holstein model of 
spinless fermions, starting from the exact atomic-limit Green function. The results agree well with 
previous calculations within the many-body coherent potential approximation. 
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I. INTRODUCTION 

Spectral properties, such as the one-electron spectral 
function, provide valuable insight into the usually com- 
plex physics of strongly correlated systems. However, 
reliable results for such quantities are difficult to obtain. 
Analytical methods are often restricted to very simple 
limiting cases, and results can usually not be extended 
to more general situations. Two remarkable exceptions 
are the Holstein model with linear electron dispersion^ 
and the Hubbard model^ both in one dimension, which 
have been solved exactly. To study more general mod- 
els, numerical methods such as Exact Diagonalization 
(ED) and Quantum Monte Carlo (QMC) have received 
much attention over the last decades. ED methods al- 
low very accurate calculations of ground-state as well as 
finite-temperature spectral properties, but are restricted 
to rather small clusters due to the large dimension of the 
corresponding Hilbert space. QMC can be used to obtain 
results on large clusters even in higher dimensions, but 
here the sign-problem and the ill-posed analytic contin- 
uation required to obtain dynamic correlation functions 
for real frequencies are detrimental for many interesting 
applications. 

The recently developed cluster perturbation theory^ 
(CPT) marks an important improvement of the situa- 
tion. It is based on a break-up of the infinite system 
into identical, finite clusters, on which the one-electron 
Green function is calculated exactly. Then, the hopping 
between clusters is treated within strong-coupling per- 
turbation theory 6,7 (SCPT). The method has been suc- 
cessfully used for various Hubbard models, to calculate 
spectral functions, as well as other quantities of inter- 
est, both for zero^*^ and finite temperature^ Although 
the concept of CPT relies on a model with local inter- 
actions only, it has also been applied with some success 
to the t — J modelii2*ii For the calculation of the cluster 
Green function, the ED Lanczos method (for a review see 
Ref.GJ 

can be used. 

For the case of coupled electron-phonon systems, such 
as the Holstein model and its various extensions-e.g. 
the Holstein-Hubbard or the Holstein double-exchange 



model- the application of ED methods is hampered by 
the infinite number of possible phonon configurations, 
which gives rise to a rapidly growing requirement of 
computer memory and/or CPU time as the number of 
lattice sites or phonon states increases. Consequently, 
standard ED methods-employing some kind of Hilbert 
space truncation-are restricted to very small clusters, es- 
pecially for low phonon frequency and/or strong electron- 
phonon coupling. Again improved methods such as 
DMRG, or the use of variational phonon bases allow 
to extend the accessible parameter range. Neverthe- 
less, as electron-phonon interaction has been identified 
as an important ingredient in, e.g., high-temperature 
superconductors^ and manganitesia, further progress 
along these lines is highly desirable. 

In this paper, we show that CPT can be successfully 
applied to electron-phonon models with a (local) cou- 
pling of the Holstein type^ We present results for the 
one-electron spectral function of the Holstein polaron, 
i.e. the Holstein model with one electron, in one and two 
dimensions. The Holstein polaron problem has been in- 
vestigated intensively in the past, and a wealth of infor- 
mation about its spectral properties is available. We find 
that the use of CPT strongly reduces finite-size effects, 
giving results which are much closer to the thermody- 
namic limit than the corresponding ED data. Addition- 
ally, we consider the special case of a completely satu- 
rated ferromagnetic state at zero temperature in the Hol- 
stein double-exchange modeU^ for colossal magnetoresis- 
tive manganites. The latter is then equivalent to the 
Holstein model of spinless fermions, and we combine the 
exact atomic-limit one-particle Green function with CPT 
for a single-site cluster to calculate the spectral function 
at quarter-filling. The results of this simple approach 
agree well with the previously developed many-body co- 
herent potential approximationii2ii£iiLi2ii2i 

The paper is organized as follows: In Sec. [H] we give a 
review of CPT. Sec. IIIII discusses the application to the 
Holstein polaron, while the SCPT for the many-electron 
case is presented in Sec. II VI Finally, Sec. IVl contains our 
conclusions. 
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II. CLUSTER PERTURBATION THEORY 



tion at zero temperature 



The basic idea of CPT^ is to divide the infinite lattice 
into identical clusters, each containing N lattice sites. 
Adopting the notation of Ref. the Hamiltonian of the 
system is written in the form H = H + V where 



G ab (e)= (n\c a 



H = J2<, V 



R 



E 

Ft, Ft' 

a.b 



-vrR.R *j* 

V a,b C Ra C R'b 



(1) 



In the most general formulation of CPT^ the subscripts 
a, b denote different orbitals within a cluster, but here we 
restrict ourselves to the case of one orbital per site so that 
a, b = 1, . . . , N. The vectors R, R' correspond to sites in 
the superlattice of clusters (see Ref. 0). In Eq. JU, Hq- 
represents a Hamiltonian of a single cluster-containing 
local interactions only-and V describes the hopping be- 
tween clusters, i.e. the hopping amplitude between site a 
of cluster R and site b of cluster R' is given by the ma- 
trix element Vf^ R . Although long-range hopping can 
also be included^ we shall only consider models with 

R R' 

nearest-neighbor hopping so that V a b ' = —t for neigh- 
boring sites a, b in adjacent clusters R, R' . Within CPT, 
an approximation for the Green function of the original 
system, (J(fc,e), is obtained using an analytical strong- 
coupling perturbation expansion up to first order of the 
inter-cluster hopping V (for details of the derivation see 
Ref. 0) . The resulting equation relating the Green func- 
tion of the original lattice to the energy-dependent cluster 
Green function G(z) reads^ 



Gab{Q,z) 



G(z) 



1 - V(Q)G(z) 



(2) 



a , b 



Here z = e + ir/, G(z) and V(Q) stand for N x N matri- 
ces, and the inter-cluster hopping V has been partially 
Fourier-transformed exploiting the translational symme- 
try of the cluster superlattice, with Q being a wave vector 
of the reduced Brillouin zone£ Finally, the Green func- 
tion Q a b can be transformed from the mixed representa- 
tion of Eq. J3J| , real space within a cluster and reciprocal 
space between clusters, using 5 . 



1 N 

Q{k,z) = - £ Gab(k,z)e- ik < r ^ 



(3) 



a,fc=l 



to obtain the familiar representation of the one-electron 
Green function. Q(k,z) as given by Eqs. © and © be- 
comes exact in the atomic limit t — (Refs. ■ More- 
over, it also reduces to the exact result for the case of 
non-interacting electrons^ since, in this case, Eq. (J2J) 
corresponds to the exact resummation of the perturba- 
tion series. Finally, CPT also becomes exact in the limit 
N — > oo (Refs.Ua). The one-electron cluster Green func- 



z - (H - E ) b 
1 



Gain) 



(4) 



can be calculated exactly for any pair of site indices a, b 
in the cluster using, e.g., the Lanczos method. Here E 
is the energy of the ground state of the cluster, and 
a spin index has been suppressed in the notation. The 
two parts of the Green function matrix G a b correspond 
to adding or removing an electron to/from |n). Finally, 
the one-electron spectral function is defined as 



A(k,e) = -% 



lim Imf/(fe,e 

?j-»0+ 



•ir?). 



(5) 



Since CPT is based on a perturbation expansion in 
the inter-cluster hopping, the method can be expected 
to work especially well in the strong-coupling regime. 
This is also illustrated by the fact that it becomes ex- 
act in the atomic limit, as mentioned above. On the 
other hand, for weak or intermediate coupling, the elec- 
tronic kinetic energy is not small compared to the local 
interactions. Consequently, the size of the cluster has 
to be large enough in order to obtain accurate results. 
In fact, from previous applications of CPT, e.g., to the 
one- and two-dimensional Hubbard model^*^ the clus- 
ter size N emerged as the main control parameter of the 
method. In the case of the one-dimensional Hubbard 
model, for example, N = 1 is identical to the Hubbard 
I approximation, 20 while N — 2 already gives a spectral 
function that contains most of the relevant features such 
as short-range antiferromagnetic ordering. 5 With increas- 
ing N, the CPT Green function approaches systemati- 
cally the exact result for the infinite system. For identical 
cluster size, the CPT spectrum contains many more poles 
with significant residues than the corresponding results 
of ED. In fact, also in the ID Hubbard model, the spec- 
trum obtained with CPT on a four-site cluster is already 
comparable in quality to the ED spectrum for N = 12 
(Refs. 00). An additional advantage of CPT is the possi- 
bility to evaluate A(k, e) at continuous wavevectors k, in 
contrast to ED which restricts k to the N vectors of the 
first Brillouin zone of the cluster, of which only N/2 + 1 
are physically distinct. Finally, finite temperature Lanc- 
zos methods can also be combined with CPT to calculate 
thermodynamic properties^ 

Concerning the application of the Lanczos method to 
calculate the cluster Green function G a b, it is important 
to stress the need for open boundary conditions (BCs). 
Attempts have been made to use periodic BCs and sub- 
tract the corresponding terms afterward in the pertur- 
bative treatment of the inter-cluster hopping, but it has 
been found that the accuracy of the results is much bet- 
ter for the case of open BCs. Although the latter are 
physically more intuitive in connection with CPT, the 
calculation of the cluster Green function with the Lanc- 
zos method becomes more difficult as one cannot exploit 
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translational symmetry. Other symmetries such as the 
inversion group can be used in principle, but are usually 
not as effective in saving computer memory by reducing 
the size of the corresponding Hamiltonian matrix and 
Lanczos vectors. 

We want to point out that CPT does not in principle 
rely on the ED method. In fact, the cluster Green func- 
tion may be calculated using any method available^ In- 
deed we will see in Sec. |W] that it is possible to combine 
the exact analytic solution for the atomic-limit Green 
function with CPT, to obtain results which agree sur- 
prisingly well with the many-body coherent potential 
approximation^ 

In addition to the spectral function considered here, 
other physical properties of the system can also be cal- 
culated with CPT. This includes, e.g., the ground state 
energy of the infinite system, the electronic kinetic energy 
or the Fermi surfaced The strength of CPT lies in the 
calculation of the one-particle Green function and related 
quantities such as the density of states. The numerical 
effort is relatively small compared to more sophisticated 
methods like DMRG or QMC. An additional advantage 
of CPT is the fact that it can easily be applied also to 
two-dimensional systems, in contrast to, e.g., DMRG. 
Finally, an important disadvantage of CPT should be 
mentioned: Within the current formulation, two-particle 
Green functions cannot be calculated. Consequently, it 
is not possible to compute, e.g., the dc conductivity or 
other interesting two-particle correlation functions. 

III. HOLSTEIN POLARON 

The Hamiltonian of the Holstein modeli^ reads 

(ij)<7 i 

-g^mibl + bt), (6) 

i 

where c\ a (c ia ) and b\ are creation (annihilation) 
operators for an electron with spin a and a phonon of 
frequency u> at lattice site i, respectively. The elec- 
tron occupation number is defined as n, = with 

riic = cJ^Cjo-, and the parameters of the model are the 
hopping integral for nearest-neighbor hopping, t, and the 
electron-phonon coupling strength, denoted as g. It is 
common to define a dimensionless coupling parameter 
A = g 2 /(2uW), where 2W is the bandwidth of the bare 
electron band, and a dimensionless phonon frequency 
Cu = uj/t. The Holstein model can then be described using 
only these two parameters. Moreover, we shall express all 
energies in units of t. As mentioned above, here we only 
consider the one-electron limit of Hamiltonian JHJl which 
is also called the Holstein polaron problem. Although 
there is only a single electron in the system, the coupling 
to the phonons makes it a complex many-body problem, 
which has been the focus of much theoretical work. The 



restriction to one electron greatly simplifies calculations 
with the Lanczos method since both, the number of re- 
quired phonon states^ and the number of electron con- 
figurations grow noticeably with the number of particles. 
However, in Sec. IIVI we will use the exact result for the 
atomic-limit Green function and CPT for a single-site 
cluster to calculate the spectral function of the Holstein 
model of spinless fermions at quarter-filling. 

Following other authorSj 22 i 23 ' 24 i 25 ' 26 i 27 we calculate the 
Green function 

S(fe,e) = (0| Cfe7 -^4|0> , (7) 

where |0) represents the ground state of the phonons and 
the vacuum state for the electrons. The spin index can 
be suppressed due to the symmetry of the problem. The 
corresponding one-electron spectral function is given by 
Eq. (5J. 

Compared to the class of Hubbard models for which 
CPT has been originally developed, we are facing an ad- 
ditional difficulty arising from the a priori infinite num- 
ber of allowed phonon states. We employ a widely used 
truncation schemed of the phonon Hilbert space which 
is spanned by the basis states 

N 1 (r) 

(r) 

where v\ denotes the number of phonons at lattice site 
i. Now the truncation consists in restricting the basis 
states to the subset with 

N 

E ^ < N ^ (») 

1=1 

leading to (N ph + N - l)!/(A ph !(A - 1)!) allowed phonon 
configurations. The convergence of the results with A p h 
can be monitored using the ground-state energy E of the 
cluster with open BCs. In all results of this paper, A p h 
was chosen such that the relative error for the ground 
state with one electron, |£;(iV p h + 1) - E(N ph )\/\E(N ph )\ , 
was smaller than 10~ 5 . We find that convergence of E 
also ensures a well-converged spectral function. More- 
over, the influence of the number of phonons kept in the 
calculation is much larger for the incoherent part of the 
spectrum than for the coherent, low-energy quasi-particle 
peak which determines E (see Sec. 1111 AT) . Finally, a re- 
fined truncation scheme which allows for extremely ac- 
curate results (relative error < 10 ) has been proposed 
by Wellein et aJA 

Before we come to a discussion of the results obtained 
with CPT, we want to comment on some of the existing 
work on spectral properties of the Holstein polaron. As 
indicated before, the most reliable method to calculate 
dynamic quantities, such as A(k, e), is ED which has been 
used extensively in the past. 22 i 23 i 24 i 25 i 26 i 27 i 28 i 29 i 30 Most of 
this work has focused on the polaron band structure E(k) 
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instead of the spectral function, since it is often easier to 
interpret, especially in the strong-coupling regime where 
the structure of A(k, e) is rather complicated. However, 
as pointed out by Wellein et al.j2& the two quantities are 
closely related. In fact the position of the lowest-energy 
peak in A(k,e), obtained from the Green function JZJ, 
follows exactly the polaron band structure as we vary 
k. Moreover, as discussed by Wellein and Fehske^S the 
integral over this first peak is equivalent to the quasi- 
particle (QP) weight z(k) = |(^|4|0}| 2 , where 
denotes the lowest-energy single-polaron state in the sec- 
tor with total momentum k. Other numerical methods 
which have been used to calculate the spectrum of the 
Holstein polaron include DMRG 3 ^ 3 ^ (in one dimension), 
finite-cluster strong coupling perturbation theory^ (ID, 
2D), QMC'''" 3 ' (1D-3D), and variational methoda 36 i 37 i 38 
(1D-4D). 



A. Comparison with Exact Diagonalization 

As mentioned above, the critical parameter of CPT is 
the number of sites in the cluster. To demonstrate the 
advantage of CPT over the standard ED method (see, 
e.g., Ref. we present in Fig. ^ the spectral function 
A(0, e) in one dimension for different cluster sizes N. We 
chose Co = 2 and A = 0.5, which is the regime where 
an extended polaron exists (see, e.g., Ref. 1231 . Conse- 
quently, significant finite-size effects can be expected for 
small clusters, which is exactly what we see in the ED re- 
sults. For the latter periodic BCs have been used. Fig. 2] 
clearly shows that the shape of the large QP peak at 
e « —2.4 changes very little with increasing N for both, 
ED and CPT, but a noticeable shift can be observed in 
the case of the ED spectra as we go from TV = 2 to N = 4. 
The influence of N is much larger for the incoherent part 
of the spectrum, which lies about a distance ui above the 
QP peak. The ED spectra display sharp, well-separated 
peaks, whereas the corresponding CPT data-containing 
many more poles-resembles much closer the expected re- 
sults for an infinite system. The latter has been investi- 
gated by Marsiglic 23 using Migdal-Eliashberg theory. For 
the same parameters, he found that the QP peak remains 
almost unchanged as N — > 00, while the incoherent part 
evolves into a continuous band that fits well to the CPT 
results even for rather small clusters N > 6. We have 
also compared A(k, e) for k ^ 0, and the observed influ- 
ence of finite-size effects agrees perfectly with previous 
work of Wellein and Fehske: 29 As k increases from k = 
to A; = it, the size of the polaron increases, and the de- 
viations of the ED data from the CPT results become 
larger. In the strong-coupling or small-polaron regime, 
not shown here, finite-size effects are known to be small. 
Consequently, even for very small clusters, ED and CPT 
both give well-converged results for the QP peak which 
determines, e.g., the ground-state energy. However, in 
the case of ED, the incoherent part of the spectrum for 
wavevector k, corresponding to excitations of an electron 
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FIG. 1: Comparison of the spectral function A(0, e) of the 
ID Holstein polaron obtained with ED (left column) and CPT 
(right column) for various numbers of lattice sites N in the 
cluster. The plot is for uo = 2.0, A = 0.5 and N p ^ = 6. An 
artificial imaginary part 77 = 0.02£ has been used to broaden 
the delta peaks. 



with momentum q and a phonon with momentum k — q, 
still exhibits the typical multi-peak structure of a finite 
system, whereas the CPT results again reproduce much 
better the incoherent band found in the thermodynamic 
limit. Moreover, as mentioned in Sec.[Hj CPT allows to 
calculate A(k, e) for continuous k, while ED on a iV-site 
cluster is restricted to N/2 + 1 physically non-equivalent 
wavevectors. 

A closer look at the CPT results in Fig.^reveals small 
additional peaks-not present in the ED spectra-which 
move from the incoherent part of A(k, e) towards the 
QP peak with increasing N. Additional calculations for 
larger clusters have shown that these peaks vanish sys- 
tematically with increasing N, so that the CPT spectrum 
approaches the exact result in the thermodynamic limit 
N = 00, as expected. Consequently, these peaks are not a 
defect of CPT, but represent finite-size effects which arise 
from the approximate treatment of inter-cluster hopping. 
The latter, in combination with the open BCs used to cal- 
culate the cluster Green function, leads to a system which 
does not have perfect translational symmetry. The situ- 
ation is equivalent to ED with open BCs: For N — » 00 
the spectrum approaches the results of an infinite cluster. 
However, in contrast to CPT, the effects for finite N are 
much more significant. Moreover, these finite-size effects 
manifest themselves in a slightly different way than in the 
case of periodic BCs, where no additional peaks-showing 
the afore mentioned behavior-are found. In the case of 
CPT, already for the small cluster sizes shown in Fig. ^ 
the spectral weight of these additional peaks is extremely 
small compared to the rest of the spectrum. For other 
values of u> and A, a similar behavior has been found. Al- 
though not discussed by the authors, similar effects can 
also be expected for the case of the Hubbard model^*^ 
although they may be larger for the Holstein polaron 
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FIG. 2: Top: Spectral function A(k, e) of the ID Holstein 
polaron calculated with CPT for N = 14, iV p h = 6 and r\ = 
0.02i. Bottom: Density plot of the same data for 100 points 
in k space. Symbols represent results of Bonca et al^S 

due to the higher sensitivity of phonon-excitations to the 
BCs. 



B. Results: One dimension 

In one dimension, the general picture emerging from 
previous work on the Holstein polaron problem is as fol- 
lows (see e.g. Ref. |28j and references therein): In the 
non-adiabatic regime (p > 1), a so-called Lang-Firsov 
polaron is formed which, due to the instant response 
of the phonons to the electronic motion, represents a 
very localized object. As the electron-phonon coupling 
increases, its mobility or effective hopping amplitude 
exhibits a gradual decrease, and for strong coupling a 
nearly-localized small polaron moving in an exponentially 
narrow band exists. In contrast, in the adiabatic (uj < 1), 
weak-coupling regime the electron drags with it an ex- 
tended cloud of phonons. This object is usually called 
a 'large polaron'. At A rs 1, a sharp but continuous 39 
transition to a less mobile small Holstein polaron takes 
place. The two conditions for small-polaron formation, 
independent of the value of uj, are A > 1 and X/uj > 0.5. 28 
Moreover, as pointed out by Capone et al.^ the forma- 
tion of small polarons is determined by A > 1 for uj < 1 



and by X/uj > 0.5 for uj > 1. Here we restrict ourselves to 
the most interesting regime of phonon energies compara- 
ble to the electronic hopping, i.e. uj ~ 1. For intermediate 
uj and A, no reliable analytical methods exist, so that nu- 
merical approaches represent the most important source 
of information. 

In Fig. |3 we present results for A(k, e) for u) = 0.8, 
A = 0.25 and N = 14, as well as a density plot of the 
same data. As mentioned before, the spectrum consists 
of a low-lying QP peak and an incoherent part at higher 
energies. The physics behind the observed behavior of 
A(k, e) has been discussed, e.g., by Stephan^S and is typ- 
ical for electronic systems weakly interacting with disper- 
sionless optical phonons. For small k, most of the spec- 
tral weight resides in the QP peak which corresponds to 
a weakly-dressed electron. For the case considered here, 
in which the phonon energy lies inside the bare electron 
band, electron and phonon hybridize and repel each other 
near the point where they would be degenerate, i.e. for 
\E(k) — E(0)\ ~ uj. This coincides with the region where 
the flattening of the polaron band occurs and, in fact, 
for larger k the phonon becomes the lowest-energy ex- 
citation. However, most of the spectral weight is con- 
tained in the broad, incoherent band which follows the 
free-electron dispersion. The density plot in Fig. [21 also 
contains data for the polaron band structure E(k) which 
has been obtained by Bonca et using their varia- 

tional ED method. The latter has been shown to give 
very accurate results for the infinite system, although it 
becomes somewhat less accurate in the strong coupling 
regime and for large values of k (Ref. .'57}. As mentioned 
before, E{k) corresponds to the lowest-energy band in 
A{k, e) and we find very good agreement with our data 
throughout the Brillouin zone. 

Fig.|3|shows results for a similar phonon frequency uj = 
1.0 but for stronger electron-phonon coupling A = 0.5 and 
N = 12. Compared to the weak-coupling case discussed 
above, the polaron band is separated more clearly from 
the incoherent part of the spectrum and, as expected, the 
band-width is further reduced. Additionally, even more 
spectral weight has been transfered to the high-energy, 
incoherent band. On top of that, a gap shows up in the 
upper band at about k — n/2. Again the polaron band 
fits very well the results for E(k) of Bonca et ali^l 

We next consider the case of intermediate coupling 
A = 1.0, with uj — 1 and N = 8 (Fig. gj. For these pa- 
rameters, an extended polaron exists which still has a rel- 
atively large band-width, compared to the small-polaron 
case discussed below. Moreover, the incoherent part of 
the spectrum has split up into several sub-bands sepa- 
rated in energy by uj, which correspond to excitations 
of an electron and one or more phonons. As before, we 
find very good agreement between the low-energy band 
in A(k, e) and the polaron band energy E{k) calculated 
by Bonca4i 

Finally, in Fig. |SJ we report the spectral function for 
uj = 1 and A = 2.0. The results have been obtained us- 
ing only a four-site cluster, which is sufficient to get very 
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FIG. 3: Top: Spectral function ^4(fc, e) of the ID Holstein 
polaron calculated with CPT for N = 12, N p ^ = 6 and r\ = 
0.02i. Bottom: Density plot of the same data for 100 points 
in k space. Symbols represent results of Bonca et alrf— 



FIG. 4: Top: Spectral function A(k, e) of the ID Holstein 
polaron calculated with CPT for N = 8, N p ^ = 9 and n = 
0.02i. Bottom: Density plot of the same data for 100 points 
in k space. Symbols represent results of Boncai— 



good agreement with Bonca's data for E(k), with only 
minor deviations at large values of k where finite-size 
effects are most pronounced, as discussed in Sec. 1111 Al 
This is a consequence of the predominantly local effects 
in the strong-coupling regime, which also manifest them- 
selves in terms of a very narrow polaron band. 



C. Results: Two dimensions 

To illustrate the applicability of CPT, we also calcu- 
lated the spectral function of the Holstein polaron on a 
2D cluster with N = 8, which has the shape of a tilted 
square. In contrast to CPT in one dimension, the choice 
for the shape of the cluster is not unique, and different 
clusters may lead to slightly different results. This possi- 
bility has been investigated for the Hubbard model, and 
the effect of the cluster shape on A(k, e) was found to 
be rather small. 5 For the case of the Holstein polaron, 
where the physics is dominated by local correlations, the 
influence of the geometry of the cluster is expected to be 
even smaller. 

As discussed, e.g., by Wellein et al.jSfi similar to ID, 
a small polaron is formed in two dimensions provided 
that A > 1 and X/ui > 0.5. While the behavior in the 



non-adiabatic regime (a) > 1) is only weakly affected by 
dimensionality, 28 42 important differences exist in the adi- 
abatic regime u> < 1: In contrast to the one-dimensional 
case, where a large polaron is formed for any A > 0, the 
electron remains quasi-free for A < 1, as indicated by 
an almost unaffected effective hopping amplitude. More- 
over, for the same value of u, the cross-over to a small 
polaron at A w 1 is much sharper in 2D than in ID. 

Here we simply aim to demonstrate the possibility of 
calculating the 2D spectral function with CPT. There- 
fore, we restrict ourselves to one set of parameters, 
namely uj = 2.0 and A = 0.945, which has also been 
treated using finite- cluster strong-coupling perturbation 
theory^. In contrast to standard SCPT based on the 
Lang-Firsov transformation, the latter has been shown 
to give reliable results also for intermediate A and to, 
which is a consequence of the inclusion of longer-ranged 
effects. 29,33 . While in the ID case the density plot of 
A(k, e) contains all 100 values of k used in CPT, in two 
dimensions we have used 400 points in k space. However 
only 60, lying along TMXr, are shown in Fig. EI 

From the above discussion, and for the parameters con- 
sidered here, we expect a rather broad polaron band. 
This is clearly confirmed by the spectral function shown 
in Fig. El and the lowest-energy band in our data resem- 
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FIG. 5: Top: Spectral function A(k, e) of the ID Holstein 
polaron calculated with CPT for N = 4, iV p h = 25 and r\ = 
0.02i. Bottom: Density plot of the same data for 100 points 
in k space. Symbols represent results of Bonca^i 
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FIG. 6: Top: Spectral function A(k, e) of the 2D Holstein 
polaron calculated with CPT for N = 8, N p ^ = 9 and r\ — 
0.02t. Bottom: Density plot of the same data (see text). 



bles closely the findings of Stephan (Fig. 2 of Ref . l33l) . 
In particular, similar to the one-dimensional case consid- 
ered in Sec. IIII Bl a flattening of the polaron band near 
(tt/2, 7r) is found which has also been noted by Wellein et 
alwS Above the polaron band, also similar to ID, there 
lie several other incoherent bands which correspond to 
multi-phonon excitations and are therefore separated in 
energy by Co. 

In summary, the results of this section clearly demon- 
strate that CPT is applicable not only in the strong- 
coupling regime, but also for weak and intermediate 
electron-phonon interaction. The quality of the result- 
ing spectra is superior to ED data for the same cluster 
size, and very good agreement has been found with the 
variational method of Bonca et alml in one dimension. 
Moreover, we have shown that CPT also allows accurate 
calculations of A(k, e) in two dimensions. 

IV. MANY-ELECTRON CASE 

In the last section, we have restricted ourselves to 
the Holstein model with one electron. Although CPT 
has been successfully applied, e.g., to the many-electron 
Hubbard model^*^^ the electron-phonon coupling in 



the Holstein model greatly complicates calculations us- 
ing Lanczos ED. For finite electron density, we combine 
CPT with the exact analytic result for the Green func- 
tion in the atomic limit. For the atomic limit, the Green 
function has been obtained for many models using the 
equation-of-motion method^ and here it will allow us 
to obtain results for the many-electron case which will 
be compared with the many-body coherent potential ap- 
proximation discussed below. 

A. Many-body coherent potential approximation 

Extending previous work of Edwards et aliZ^ for the 
pure double-exchange (DE) model (see, e.g., Ref. 0), 
Greeniii studied the Holstein-DE model using a many- 
body coherent potential approximation (CPA) which, 
owing to the more complicated form of the Holstein-DE 
Hamiltonian, constitutes a considerable extension of the 
Hubbard III approximation^ The many-body CPA suc- 
cessfully describes many aspects of the manganites, and 
we refer the reader to a recent review of this work by 
Edwards ^ Here we only consider the special case of a 
completely saturated ferromagnetic state at temperature 
T — 0, with all itinerant spins having f spin, say. Con- 
sequently, the DE term which couples local and itinerant 
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spins 1 * becomes merely a constant shift in energy, and 
the Holstcin-DE model is equivalent to the pure Holstcin 
model of spinless fermions, i.e. with no doubly-occupied 
sites' An important feature of the many-body CPA is 
that the one-electron Green function reduces to the exact 
atomic limit for t = 0, which takes the formic 

Gf L ( e ) = e-4i + V^P^ + !^)), 
1 w It t-^f r\ \e + Lor e-ur) J 

(10) 

where a = g 2 /uj 2 and the polaron binding energy 
— (g 2 /cu)n (n — 0, 1) has been absorbed into the chemi- 
cal potential. The general result for G AL of the Holstein 
model with electrons of both spins has been given by 
Green^i and we drop the spin-index in the sequel. As 
discussed by Edwards?^ for an elliptic density of states, 
the local Green function Q(z) for complex energy z sat- 
isfies the CPA equation 

g{z) = G AL (z~w 2 g/4) (ii) 

and the self-energy can be obtained from 1 ^ 

e(z) = z-g- 1 -w 2 g/4. (12) 

Finally, the one-electron spectral function is given by 

A{k,e) = -TT- 1 lm[z-e k -^{z)] 1 (13) 

where 

3 

e k = -2t cosfc m (14) 

m— 1 

is the band energy for wavevector k. 

In order to compare with angle-resolved photoe- 
mission (ARPES) data on the bilayer manganite 
ha 1 2 Sr 1 8 Mn 2 7 , nominally with n = 0.6, Hohenadler 
and Edwards chose a strong-electron phonon coupling 
g/W — 0.2, as deduced from the low Curie temperature 
of this material^ To simplify calculations, they also used 
n = 0.5 for which case the chemical potential /i = by 
symmetry. We want to point out that the many-body 
CPA assumes a homogeneous system, so that no ten- 
dencies toward charge-density-wave order occur as n is 
varied." As in previous workjiSi Hohenadler and Edwards 
used W = leV and u/W = 0.05 (see also Ref. El). The 
results-^ for A(k, e), shown in Fig. d support the theory 
of Alexandrov and Bratkovsk}^ that in these mangan- 
ites, small polarons exist in the ferromagnetic state. A 
similar interpretation of the experimental data-based on 
standard small-polaron theory-had also been given by 
Dessau et ali£ Well away from the Fermi surface, a well- 
defined peak exists which broadens as k approaches the 
Fermi level E F at y = 0.5. If y is increased further, most 
of the spectral weight is transfered above E F . More- 
over, the peaks never approach the Fermi level closely, 
in agreement with the experimental data. This indicates 
the existence of a pseudogap in the one-electron density 
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FIG. 7: Comparison of the spectral function of the Holstein 
model of spinless fermions at T = 0, calculated with SCPT 
(top ) an d with the many-body CPA (bottom, taken from 
Ref. Il9t) . Here the wavevector k is given by k = 7r(l,y,0) 
with y as indicated in the figure. The plot is for u/W = 0.05 
and g/W — 0.2. The SCPT results have been broadened 
using a smearing parameter r//W = 0.05. 

of states. However, in the gap, there exist small polaron 
sub-bands (see Fig. 4 of Ref. fig) and one of them, at the 
Fermi level, presumably gives rise to the low but finite 
conductivity of the system. As discussed by Edwards, 14 ^ 
the many-body CPA does not give coherent states with 
infinite lifetime at the Fermi level, even for T = 0. This 
is typical for any CPA, and here it leads to an incoherent 
polaron sub-band around the Fermi level. Nevertheless, 
outside the central band around E F , the imaginary part 
of the self-energy displays the correct behavior, i.e. it 
vanishes in the gap, between the polaron bands. 

B. SCPT 

In this section we use the exact result for the atomic- 
limit Green function of the Holstcin model of spinless 
fermions, G AL (Eq. {TO)), and combine it with CPT 
to compare the resulting spectrum with the many-body 
CPA. For this case of a single-site cluster (a = b), Eqs. @ 
and JSJ reduce to a single equation for the one-electron 
Green function^ 
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with z = e + \r) and e k as defined by Eq. (|14fl . Hence, 
as mentioned before, CPT for N = 1 is equivalent to 
the Hubbard I approximation^ but here with the more 
complicated atomic-limit Green function of the Holstcin 
model given by Eq. 1|10|) . In the sequel, we shall re- 
fer to this approximation as SCPT. This is justified by 
the fact that the approach becomes exact for t = 0. 
Historically a similar strong-coupling expansion for the 
Hubbard model£i£-including higher order corrections- 
has been the starting point for the development of CPT. 

Before we discuss the results, we would like to comment 
on the quality of the SCPT used here: While the many- 
body CPA requires a self-consistent, iterative solution 
of Eq. (U , the SCPT Green function is obtained from 
the Lehmann representation of the atomic-limit Green 
function ijTUI) . and the subsequent use of the resulting 
self-energy S AL in Eq. (|15fl . Similar to the original Hub- 
bard I approximation^ the resulting Green function con- 
sists of delta peaks corresponding to states with infinite 
lifetime. However, due to the poles in the self-energy, 
there are no states at the Fermi level and the system is 
not a Fermi liquid. As in the many-body CPA, Q de- 
pends on k only through the band energy e fc , whereas 
the self-energy is local. This reliance on the atomic limit 
is reasonable in the strong-coupling regime considered 
here, where small polarons move in an extremely narrow 
band. Consequently, the simple perturbative treatment 
of the hopping term can be expected to give sensible re- 
sults. Nevertheless, in SCPT, we have to use an artificial 
imaginary part 77-which does not depend on energy-to 
obtain peaks of finite width. Although for large enough 
77 there will be states at the Fermi level, the latter have 
only finite lifetime even for T — 0. Hence, both SCPT 
and the many-body CPA never give a Fermi liquid, but 
the self-consistent CPA Green function yields an imagi- 
nary part of the self-energy that shows the correct, strong 
energy-dependence except for the region inside the very 
small, incoherent polaron band around E Fl as discussed 
in Sec. IIV Al Thus, as could be expected from the Hub- 
bard I-like approximation in Eq. (JT3J, the many-body 
CPA is superior to SCPT, although both approaches be- 
come exact in the atomic limit. 

The spectral function obtained with SCPT using 
Eq. © , also shown in Fig. resembles quite closely the 
results of Hohenadler and Edwards. 19 Although there are 
some differences concerning the width and the position 
of the peaks, the overall behavior is very similar. In par- 
ticular, the broadening of the QP peak near the Fermi 
surface at y = 0.5 is well reproduced. Clearly the success 
of SCPT consists in a surprisingly good agreement with 
the CPA data for all k. Despite this agreement, CPT 
fails to reproduce the polaron sub-bands, and the sharp 
edge to the pseudogap for large values of y. Moreover, 
the gap is larger than in the CPA data. These short- 
comings are a consequence of the rather crude approxi- 
mation. Nevertheless, keeping in mind the simplicity of 
the ansatz, the agreement with the many-body CPA is 
satisfactory We would like to point out that the SCPT 



presented here can also be generalized to the Holstein-DE 
model with quantum spins (e.g. 5 = 3/2 appropriate for 
the manganites 1 ^) and at finite temperature, using the 
atomic- limit Green function given by Edwards. 14 Finally, 
the approximation could be systematically improved by 
increasing the number of sites in the cluster, which is 
exactly the idea behind CPT. 

However, for N > 1 the cluster Green function can no 
longer be calculated analytically and one has to resort to 
numerical methods such as ED as in Sec. IIIII Such cal- 
culations are extremely difficult for the case of quarter- 
filled two- or three dimensional clusters, small phonon 
frequency and strong electron-phonon coupling. Fu- 
ture work along these lines-employing optimized phonon 
approaches 4 ^ (see also Sec.EJ-is highly desirable in order 
to assess the quality of the many-body CPA results. 



V. CONCLUSIONS 

We have applied cluster perturbation theory to the 
Holstein polaron problem in one and two dimensions, and 
comparison with existing work has revealed very good 
agreement. In combination with the Lanczos method to 
calculate the cluster Green function, the method gives 
reliable results for the one-electron spectral function 
A(k,e) 1 which become exact in the weak- and strong- 
coupling limit, A = and t = 0, respectively, and for the 
case of an infinite cluster. Calculations for continuous 
values of the wavevector k are possible and, more impor- 
tantly finite-size effects are significantly reduced com- 
pared to standard ED. Our results extend previous ap- 
plications of CPT to Hubbard and t — J models, showing 
that the method is also well-suited for electron-phonon 
models with local interactions. As pointed out before, 
for more than one electron in the system, it becomes in- 
creasingly difficult to include enough phonon states so as 
to obtain converged results. Future work may therefore 
combine optimized phonon approaches 47 and CPT to in- 
vestigate more complicated problems such as the many- 
electron case, or extended models including, e.g., a Hub- 
bard term. The major advantage of such an approach 
is the possibility of using relatively small clusters while 
still obtaining results which are only weakly influenced 
by finite-size effects. 

Additionally we have used the exact atomic-limit 
Green function of the Holstcin model of spinless fermions, 
to calculate the spectral function for the case quarter- 
filling and strong electron-phonon coupling. The results 
of this approximation at the Hubbard I level agree sur- 
prisingly well with the many-body CPA. 
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